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, <^ ■ Abstract 

In this paper wc investigate a class of swing options with firm constraints in view of the 
' modehng of supply agreements. We show, for a fully general payoff process, that the premium, 

solution to a stochastic control problem, is concave and piecewise affine as a function of the 
P^ '. global constraints of the contract. The existence of bang-bang optimal controls is established 



X 



for a set of constraints which generates by affinity the whole premium function. 

When the payoff process is driven by an underlying Markov process, we propose a quanti- 
zation based recursive backward procedure to price these contracts. A priori error bounds are 



. established, uniformly with respect to the global constraints. 
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^ ■ 1 Introduction 

C5 I The deregulation of energy markets has given rise to various famihes of contracts. Many of them 
' appear as some derivative products whose underlying is some tradable futures (day-ahead, etc) on 
gas or electricity (see [12] for an introduction). The class of swing options has been paid a special 
I attention in the literature, because it includes many of these derivative products. A common 
feature to all these options is that they introduce some risk sharing between a producer and a 
trader, of gas or electricity for example. From a probabilistic viewpoint, they appear as some 
stochastic control problems modeling multiple optimal stopping problems (the control variable 
I is the purchased quantity of energy); see e.g. [lOl [9] in a continuous time setting. Gas storage 
contracts (see [6], [8]) or electricity supply agreements (see [18], |7]) are examples of such swing 
options. Indeed, energy supply contracts are one simple and important example of such swing 
options that will be deeply investigated in this paper (see below, see also |12j for an introduction). 
It is worth mentioning that this kind of contracts are slightly different from multiple exercises 
American options as considered in [10] for example. In our setting the volumetric constraints play 
a key role and thus, the flexibility is not restricted to time decisions, but also has to take into 
account volumes management. 

Designing efficient numerical procedures for the pricing of swing option contracts remains a 
very challenging question as can be expected from a possibly multi-dimensional stochastic control 
problem subject to various constraints (due to the physical properties of the assets like in storage 
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contracts). Most recent approaches developed in mathematical finance, especially for the pricing 
of American options, have been adapted and transposed to the swing framework: tree (or "forest" ) 
algorithms in the pioneering work [T7], Least squares regression MC methods (see [B]), PDE's 
numerical methods (finite elements, see |29|). 

The aim of this paper is to deeply investigate an old question, namely to elucidate the structure 
of the optimal control in supply contracts (with firm constraints) and how it impacts the numerical 
methods of pricing. We will provide in a quite general (and abstract) setting some "natural" 
(and simple) conditions involving the local and global purchased volume constraints to ensure the 
existence of bang-bang optimal strategy (such controls usually do not exist). It is possible to design 
a priori the contract so that their parameters satisfy these conditions. To our knowledge very few 
theoretical results have been established so far on this problem (see however [6] in a Markovian 
framework for contracts with penalized constraints and [27j . also in a Markovian framework). 

This first result of the paper not only enlightens the understanding of the management of a 
swing contract: it also has some deep repercussions on the numerical methods to price it. As a 
matter of fact, taking advantage of the existence of a bang-bang optimal strategy, we propose and 
analyze in details (when the underlying asset has a Markovian dynamics) a quantized Dynamic 
Programming procedure to price any swing options whose volume constraints satisfy the "bang- 
bang" assumption. Furthermore some a priori error bounds are established. This procedure turns 
out to be dramatically efficient, as emphasized in the companion paper [5] where the method is 
extensively tested with assets having multi-factor Gaussian underlying dynamics and compared to 
the least squares regression method. 

The abstract swing contract with firm constraints The holder of a supply contract has the 
right to purchase periodically (daily, monthly, etc) an amount of energy at a given unitary price. 
This amount of energy is subject to some lower and upper "local" constraints. The total amount 
of energy purchased at the end of the contract is also subject to a "global" constraint. Given 
dynamics on the energy price process, the problem is to evaluate the price of such a contract, at 
time t = when it is emitted and during its whole life up to its maturity. 

To be precise, the owner of the contract is allowed to purchase at times tj, i = 0, . . . , n — 1 a 
quantity qi of energy at a unitary strike price Ki := K{ti). At every date ti, the purchased quantity 
qi is subject to the firm "local" constraint, 

Qmin ^ Qi ^ Qma,xi i — 0, . . . , 7T. 1, 

whereas the global purchased quantity qn ■= "^^Zo Qi is subject to the (firm) global constraint 

] (0 < Qmin < Q max < +oo). 

The strike price process (-ftri)o<i<n-i can be either deterministic (even constant) or stochastic, 
e.g. indexed on past values of other commodities (oil, etc). Usually, on energy markets the price is 
known through future contracts (-Fs,t)o<s<t where Fg^t denotes the price at time s of the forward 
contract delivered at maturity t. The available data at time are (-Fb,t)o<t<T (in real markets this 
is of course not a continuum). 

The underlying asset price process, temporarily denoted {StJo<i<n-i, is often the so-called 
"day-ahead" contract i*t,t+i which is a tradable instrument or the spot price Ft^t which is not. All 
the decisions about the contract need to be adapted to the filtration of {St J i.e. Ti := a{Stj,j = 
0, . . . , z), i = 0, . . . , n — 1 (with JFq = {0, f^}). This means that the price of such a contract is given 
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at any time tk, by 




n-l 



Qj • (^) -^j) ^ [Qmiiii Qma,x\ij — k, ... ,71 1, ^ ^ Qj € [Qmim Qmax] 

j=k 



where Q'^-^ = Q 

min Qk: Qma.x — Qmax Qk denote the residual global constraints and r denotes the 
(deterministic) interest rate. This pricing problem clearly appears as a stochastic control problem. 

In the pioneering work by [TTj, this type of contract was computed by using some forests of 
(multinomial) trees. A natural variant, at least for numerical purpose, is to consider a penalized 
version of this stochastic control. Thus, in [6], a penalization QniVn,<ln) with Qn{x,q) = —{{x — 
Qmax)"^ + {x - Qmin)~^)/£ is added {Qn is negative outside [Qmin, Qmax] and zero inside). 

As concerns more sophisticated contracts (like storages), the holder of the contract receive a 
quantity ^{ti, Sj. , qi) when deciding Qi. When dealing with gas this is due to the storing constraints 
since injecting or withdrawing gas from its storing units induce fixed costs (and physical constraints 
(pressure, etc)). 

As concerns the underlying asset dynamics, it is commonly shared in finance to assume that the 
traded asset has a Markovian dynamics (or is a component of a Markov process like with stochastic 
volatility models). The dynamics of physical assets for many reasons (some of them simply coming 
from history) are often modeled using some more deeply non-Markovian models like long memory 
processes, etc. 

All these specific features of energy derivatives suggest to tackle the above pricing problem in a 
rather general framework, trying to avoid as long as possible to call upon Markov properties. This 
is what we do in the first part of the paper where the general setting of a swing option defined 
by an abstract sequence of ^^-adapted payoffs is deeply investigated as a function of its global 
constraints (Qmim Qmax) (when the local constraints are normalized i.e. Qi is [0, l]-valued for every 
i € {0, . . . , n — 1}). We show that this premium is a concave, piecewise afRne, function of the global 
constraints, affine on triangles of the (m, M) + {(n, v), < u < v < 1} , m, M £ N^, m < M < n 
and {m,M) + {{u,v), 0<v<u<l}, m, Me'N'^,m<M-l<n-l. We also show that 
for integral valued global constraints, the optimal controls are always bang-bang i.e. the a priori 
[0, l]-valued optimal purchased quantities q* are in fact always equal to or 1. Such a result can 
be extended in some way to any couple of global constraints when all the payoffs are nonnegative. 

Then, when there is an underlying Markov "structure process", we propose an optimal quan- 
tization based on numerical approach to price efficiently swing options. This Markov "structure 
process" can be the underlying traded asset itself or a higher dimensional hidden Markov process: 
such a framework comes out in case of multi- factor processes having some long- memory properties. 

Optimal Quantization was first introduced as a numerical method to solve nonlinear problem 
arising in Mathematical Finance in a series of papers [UElISlll] devoted to the pricing and hedging of 
American style multi-asset options. It has also been applied to stochastic control problem, namely 
portfolio optimization in [24]. The purely numerical aspects induced by optimal quantization, with 
a special emphasis on the Gaussian distribution, have been investigated in [26]. See |25j for a survey 
on numerical application of optimal quantization to Finance. For other applications (to automatic 
classification, clustering, etc), see also [H]. In this paper, we propose a quantized backward dynamic 
programming to approximate the premium of a swing contract. We analyze the rate of convergence 
of this algorithm and provide some a priori error bounds in terms of quantization errors. 

We illustrate the method by computing the whole graph of the premium viewed as a function of 
the global constraints, combining the affine property of the premium and the quantized algorithm 
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in "toy model" : the future prices of gas are modeled by a two factor Gaussian model. An extensive 
study of the pricing method by optimal quantization is carried out from both a theoretical and 
numerical point of view in [3] . 

The paper is organized as follows. In the section below we detail the decomposition of swing 
options into a swap contract and a normalized swing option. In Section O we precisely describe 
our abstract setting for normalized swing options with firm constraints and the variable of interest 
(global constraints, local constraints, etc). In Section [H we establish the dynamic programming 
formula satisfied in full generality by the premium as a function of the global constraints (this 
unifies the similar results obtained in Markov settings, see [T7], [6], etc) and we show this is a 
concave function with respect to the global constraints. Then, in Section [5l we prove in our 
abstract framework that the premium function is piecewise affine and that the optimal purchased 
quantities satisfy a "0-1" or bang-bang principle (Theorem [2]). A special attention is paid to the 
2-period model which provides an intuitive interpretation of the results. In Section [6l after some 
short background on quantization and its optimization, we propose a quantization based backward 
dynamic programming formula as a numerical method to solve the swing pricing problem. Then 
we provide some error bounds for the procedure depending on the quantization error induced by 
the quantization of the Markov structure process. 

Notations. • The Lipschitz coefficient of a function / : M'^ ^ M is defined by [/]Lip := 
sup fiy)\ ^ -j-oo. The coefficient [/]Lip is finite if and only if / Lipschitz continuous. 

1-^ y\ 

• The canonical Euclidean norm on will be denoted I . 1. 



2 Canonical decomposition, normalized swing option 

As a first step we need to normalize this contract to reduce some useless technical aspects. In 
practice this normalization, in fact decomposition, corresponds to the splitting of the contract 
into a swap and a normalized swing. The decomposition can be derived from the fact that an 
^-measurable random variable q is [(^min, 9max]-valued if and only if there exists a [0, l]-valued A- 
measurable random variable q' such that q = qmm+iQmax—QminjQ' ■ Then, for every /c G {0, . . . , n— 1}, 

n-l 

PkiQLn, Qmax) = Q^nin Y^e""^^'^ "''^^(St, " K, \ F^) HQ^.. " (?min)Pf (qLo Qmax 



j=k 



swap contract 

where 



normalized contract 



P^k _ Qmin ^)gmin , _ Qmax k)'ln 

Vmin ~ „ „ dlia Vmax 



9max 'Zmin Qmax Qi 



mm 



and p^^'^l'" ^Q^in,Qmaxj is a normalized swing contract in which the local constraints are [0, 1]- 
valued . 



3 An abstract model for swing options with firm constraints 

One considers a sequence (Vfc)o<A:<n-i of integrable random variables defined on a probability 
space {Q, A, P). Let := a{Vo, Vi, . . . , Vk), k = 0, . . . ,n — 1 denote its natural filtration. For 
convenience we introduce a more general discrete time filtration := {J'k)o<k<n-i to which V = 
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{Vk)o<k<n-i is adapted i.e. satisfying C Tk, k = 0, ...,n — 1 i.e. such that the sequence 
(VA:)o<fc<n-i is J'-adapted. 

We aim to solve the following abstract stochastic control problem (with maturity n) 

(5)^ = esssup I E I ^gfcVfc \J='o\,qk:i^,Tk)^[0,l],0<k<...,n-l,^qk£ [Qmm, Qmax] > 

I \fc=0 / it=0 J 

(3.1) 

where Qmin and Qmax are two non-negative ^o-™easurable random variables satisfying 

< Qmin < Qmax < n. (3.2) 

(The inequality Qmax < ^ induces no loss of generality: one can always replace Qmax by QmaxAn 
in (jS.ip since in any case Qo + ■ • • + Qn ^ n-)- 

Note that no assumption is made on the dynamics of the state process V = (Vfc)o<fc<n-i- 
We need to introduce the following notations and terminology: 

- Throughout the paper, esssup will always be taken with respect to the probability P so P will 
be dropped from now on. 

- A couple Q := (Qmim Qmax) of non-negative .Ffc-measurable random variables satisfying < 
Qmin ^ Qmax ^ n—k is Called a couple of global constraints at time k. 

- An .F-adapted sequence q = ((?fc)o<fc<n-i of [0, l]-valued r.v. is called a locally admissible 
control. For any locally admissible control, one defines the cumulative purchase process by 

go := 0, qf,:=qQ-\ ^ q^-i, k = 1, . . . ,n. 

If 9n£ [Qmiin Qmax] , Q is Called an {T,Q)- admissible control. 

- For every k£ {0, . . . , n} and every couple of global constraints Q := (Qmin, Qmax) at time k, 

set 

PkiQ, {V, T)) := esssup |e 1^^^ I ' ^ (•^' Q)-admissible control 

so that Pq{Q, {V,J-)) is the value function of the stochastic control problem (j3.1|) when the global 
constraints Qmin and Qmax (at time 0) satisfy (13. 2p . Note that the standard convention esssup(0) = 
yields = 0. To alleviate notations, the payoff process V the filtration will be often dropped 
inP,-(Q,(y,^)). 

To be precise we will answer the following questions: 

• Existence of an optimal control q* = {ql)o<k<n-i- 

• Regularity of the value function Q i— > PjliQ, (V, JT)). 

• Existence of a bang-bang optimal control {ql.)o<k<n-i for certain values of the global constraints 
Q (namely when Q has integral components) ? 

By bang-bang we mean that W{duj)-a.s. 

- all the local constraints on the qk{^) are saturated i.e. for every k G {0, . . . ,n — 1}, qk{uj) S 

{0, 1 A Qmax} 

or 

- there exists at most one instant ko{uj) such that qko(uj)i^) ^ (0,1'^ Qmax) and one global 
constraint is saturated. 

Note that if Q G N^, then a bang-bang Q-admissible control necessary satisfies ¥{duj)-a.s. 
qk{'^) G {0, 1}, A; = 0, . . . , n — 1. The existence of bang-bang optimal controls combined with the 
piecewise affinity of P^{Q) will be the key in the design of a numerical. 



(3.3) 
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• When there is an underlying structure Markov process {Vk = ffc(lfc)), we will show that the 
optimal control turns out to be a function of Yfc at every time k as well. 

4 Abstract dynamical programming principle and first properties 
4.1 Basic properties 

As a first step, we need to establish the following easy properties of PJ^ as a function of the global 

constraints Q = (Qmin, Qmax)- 

PI. For every k£ {0, . . . ,n}, and every couple of global constraints (at time k), 

PkiQ^ (yk, J^k)o<k<n-l) = Pq ^ iQ ^ iVk+t^ ^k+l)Q<l<n-k-l) ■ 

This is obvious from p.3p . 

P2. If Q and Q' are two admissible global constraints (at time k) then 

P^{Q) = P]^{Q') on the event {Q = Q'}. 

Proof. Owing to PI one may assume without loss of generality that /c = 0. Let q and q' be 
two admissible controls with respect to Q and Q' respectively. Set qi = g^ljQ^Q/} + q'fli^Q^Qiy, 
£ = 0, . . . , n — 1. The control q is admissible with respect to Q' since {Q = Q'} G J-q. Furthermore, 

(n-l \ /n-l \ /n-l 

^qiVi>\To\ = E ^ l|Q=Q/}g^V£ I -^^0 = 1{Q=Q'}IE ^ g^V^ | JTq 

< l|Q = Q,|Po"(Q') ci.s. 

Hence, 

1{q=q'}Pq{Q) = esssup |l{Q=Q/}E i^^qtVi> \ , q Q-admissible| < 1^q=q>}Pq{Q') a.s. 
The equality follows by symmetry. ^ 

P3. Let k(z {0, . . . , n — 1}. The set of admissible global constraints Q = {Qmin, Qmax) «^ time k is 
convex and the mapping Q PJ}{Q) is concave in the following sense: if Q and Q' are two couples 
of admissible constraints, then for every random variable A : (O, J-^) ~^ [0, 1], AQ + (1 — X)Q' is an 
admissible couple of constraints and 

P^iXQ + (1 - X)Q') > XP^iQ) + (1 - X)P^{Q') a.s. 

Furthermore Qmin ^ -Pfc((3mm,Q max) is non-increasing and Qmax ' — ^ Pk{Qmin:Qma,x) is non- 
decreasing i.e. 

P^{Q)<P^{Q') a.s. on {Q^^i^ < Q„,in < Qmax < Q^axl- 

Proof. One may assume by PI that A; = 0. The convexity of admissible global constraints is 
obvious. As concerns the concavity of the value function, note that if q and q' are locally admissible 
controls then Xq + {1 — X)q' := {Xqk + (1 — A)g^)o<fc<n-i is still locally admissible. If q and q' satisfy 



6 



the global constraints induced by Q and Q' respectively, then \q+ (1 — X)q' always satisfies that 
induced by XQ + (1 — X)Q' . Consequently, using that A is J^o-™easurable, 

Po"(AQ + (l-A)Q') 

( n-l 



> esssup I E(^(Ag'fe + (1 - >^)<l'kWk \ ^o), q, q' locally admissible, g„ G [Qmin, Qmax], q'n 

[ k=0 

= A esssup <^ IE(^ qkVk I -^o), q locally admissible, qn^ [Qmin, Qmax] > 



fc=0 

r n-l 



+ (1 — A) esssup < IE(^ q'k^kl^o), q' locally admissible, q'^ G [Qmmi Qm 
I k=0 

XP^iQ) + (1 - A)Po"(Q')- 



The monotony property is as follows: let A := {Q^i^ < Qmin < Qmax < Qmaxl ^ -^o and let g' 
be a fixed Q'-admissible control. Then, for every Q-admissible control q, set 

q' ■= q'^A + ^''^"A- 

Then q' is clearly Q'-admissible and 

\fc=0 / \fc=0 / 

/ n-l \ 



A:=0 



< l^PSm a.s. 

so that ^APoiQ) < ^APoiQ') a.s. ^ 



P4. Let k G {0, 



,n 



1}. Let Q^'^\ n > 1, be a sequence of admissible global constraints such 



that Qmax ^0 as n oo. Then 

P^"(g{«))^0 



as 



n ^ oo. 



Proof. Owing to PI, one may assume k = without loss of generality. The result is straightfor- 
ward once noticed that for every Q-admissible control 



/n-l 



E[Y,qkVk\To 



\k=0 



< QmaxE max |Vfc||J^o 







P5. Let T~^{n) := {{u,v) G M^, < u < t; < n} and let k & {0, . . . ,n — 1}. There is process 
: (T+(n) X n, fior(T+(n)) ® J^fc) ^ M such that 

(i) {u,v) I— > Il^[u,v,uj) is concave and continuous on T~^{n) for every wG 

(ii) u 1-^ Il'^{u,v,uj) is non-increasing on [0,v] and v Il'^{u,v,uj) is nondecreasing on [u,n] for 
every loQ Q. 

(Hi) For every admissible constraint Q = (Qmim Qmax) o-t time k, PJ]'{Q){u!) = II'^{Q{uj),uj) ¥{duj)- 
a.s. 
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Proof. Classical consequence of P3: for every (r,s)G r+(n) set 

UUir,s),u;) :=P^{ir,s)){Lo). 

Then, set for every {u,v) G T'^{n) 

UUiu,v),u;):=a.s lim P,"((r,.))(a;). 

{r,s)^{u,v),{r,s)eT+ {n)nQ^ ,r<u<v<s 

One shows using the concavity and monotony properties established in P3 that the above 
limit does exist and that , W{duj)-a.s. {u,v) i— > IT^{{u,v),uj) is continuous on T^{n) and that 
P,"(Q)H = n^(QH,a;)P-a.5.. 

P6. Let Q = (Qmini Qmax) be a couple of global constraints (at time 0). 
Pq{Q) = esssup < E Qe{y^i)^ \J^o\,Q i^, (5)-admissible > on the event {Qmin = 0} € J^q. 



■,£=0 



Proof. This follows from the simple remark that one can define from any (T, (5)-admissible control 
q a new (J^, (5)-admissible control qhy qi := qe'i- {Ve>o}n{Q min=o} + Qi'^{Q n.in>o} and 

n—1 n—1 n—1 

l{gmin = 0} Yl ^^^^ - l{Qmin = 0} ^ ^^(V^)^ = l{Q^i„ = 0} ^ ^^(Vf)"^- 

A;=0 A:=0 fc=0 

4.2 Dynamic programming principle 

The main consequence of P5 is that at every time A;, one may assume without loss of generality 
that the couple of admissible global constraints Q = (QmimQmax) is deterministic since, for any 
possibly random admissible constraint Q (at time k) and every 0, P^{Q){uj) = Wl{uj,Q{ijj)). 

As a consequence, for notational convenience, we will still denote instead of IIJJ so that for 
any admissible global constraints Q at time k 

P,"(Q)(^) = P,"(a;,x)|,=Q(,). 

Theorem 1 (Backward Dynamic Programming Principle) Set P^ = 0. 

(a) Local Dynamic programming formula. For every k€ {0, . . . ,n — 1} and every couple 
Q = (Qmiii) Qmax) of deterministic admissible global constraints at time k 

P^iQ) = sup [xV, + E (p,\i(x"-^-^(g,x)) I ^fc) , /^"i-'^} (4.1) 

where x^iQ,x) = ((Qmin - 2;) + , (Qmax - x) AAf) andl^^ := [(Qmin - A^)+ Al, Q„,axAl]. 

(6) Global Dynamic programming formula. For every couple Q = (Qmin, Qmax) of admissible 
global constraints at time 0, the price of the contract at time k G {0, . . . , n — 1} is given by PjJ^Q^'*) 
where 

Q^'* := {Q raax — <lkiQ rain — <lk) (residual global Constraints) (4-2) 

with ql = ql{Q^'*) 

qliQ) := argmax^g,„-i-. (xFfc+E(p,!Vi(x"-'-'(Q,^))l-^fc)), ^ = 0,...,n-l. (4.3) 

Furthermore, 

/n-l \ 
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Remark. The definition (j4.3p may be ambiguous when argmax is not reduced to a single point. 
Then, one considers min argmax to define q\{Q)- 

Proof, (a) It is clear that owing to PI the case A; = is the only one to be proved. As a first step, 
we prove that 

P^{Q) = esssup [qoVo + E (Pr(x"-'(Q, ^o)) I ^o) , go : (f^, ^o) ^ /q"'} (4.4) 
: Let q = ((?fc)o<fc<n-i be a Q-admissible control. Then, qo is [0, l]-valued (as well as the g^'s) 



< 



and 



Qmin - (n - 1) < Qmin " (^l H h qn-l) < qo < Qmax " (^1 H h Q'n-l) < Qr 



so that go is Iq ^-valued. Furthermore, 

max(0, Qniin - go) < <?! H \- q-n-l < (Qmax " ^o) A(n - 1). (4.5) 

note that x'^'^iQ, 

qo) is an admissible couple of (JFg-measurable) constraints at time 1. Conse- 
quently, 

E Yl 1i^^\^o] = E E Yl QeVelTi] \To 

yi<^<n-i J \ yi<^<7i-l J 

< E(P['(x"~'(Q,go))|-^o) a.s. 

where the last inequality follows from the definition of P", (|4.5p and the monotony of con- 

ditional expectation. Then, 



E I Yl <liVi\To\<E {qoVo + Pi"(x"^'(Q, qo)) I ^o) a.s. 
One concludes that 



^o<e<n-i 



P^iQ) < esssup {e {qoVo + A"(x""'(Q, ^o)) I ^o) , qo : {^,^"0) ^ Iq''} ■ 



> 



: We proceed as usual by proving a bifurcation property for the controls. Let qo and q^ be two 
1q ^-valued .Fo-measurable random variables. Set 

^0 := {qoVo + E (Pi"(x""'(Q, qo)) I ^0) > ^o^o + E (Pr(x"-'(Q, ^o)) I -^0) } G -^o- 

and 

qo = golAo + Qo^'^Ao ■ {^,^0) Iq~'- 

Then, one checks that 

goVb + E(A"(x"-^(g,go))|^o) = max^ (y^o + E (Pf (x"-^(g, y))|./^o)) • 

Consequently there exists a sequence q^^^ of [0, Qmax A 1] -valued random variables such that 
esssup [qoVo + E (Pr(x"-'(Q, ^o)) I -^0) , qo : -^o) ^ Iq"'} 

= sup g("Vo + E (pr(x"-Hg, gi"^)) I -^0 

n 



lim ^r^o+E pr(x"-^(Q,gr))i-^o 
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One may assume by applying the above bifurcation property that the above supremum holds as a 
nondecr easing limit as n ^ oo. 

Now, for every fixed n > 1, there exists a sequence of [0, l]"~^-valued random vectors ((7^"''"^)i<fc<n-i 

such that YJk=l Qk'''^^ ^ [(Qrnin - fZo"^)^' Qmax " ^o"^] ^"^^ 



\k=l 



where we used that the admissible sequences {qk)i<k<n~i ^or the problem starting at 1 clearly 
satisfy the bifurcation principle due to the homogeneity of conditional expectation E( . \J-'o) with 
respect to Tq -measurable r.v.. Consequently, 
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(n) 



Vo+Pnx''-\Q,qt'')) = lim\("Vo + EK;g("'™Vfc|^i a.s. 



\k=l 
n-1 



EU"Vo + Pf(x"-^(Q,g(")))|^o) = lim'EU("Vo + 5^4"'™V,|^o a.s. 



k=l 



where we used the conditional Beppo Levi Theorem. Note that for every n, m > 1, {qQ^\ q^'^\ k 
1, . . . , n - 1) is an admissible control with respect to Q since x + (Qmin — x) ^ Qm'm- Hence 



n-1 



E (gi"Vo + Pi"(x""H'3, 'zS"^)) I -^o) < esssup |e \^qoVo + Y^mVk | j , <? Q-admissible| a.s. 

= P^iQ) a.s. 

To pass from (Of) to (liTTD is standard using P5. Let qo : {0,,To) Iq~^. 

(?o^o + E(Pr(x"-^(Q,go))|^o) = go^^o + E(Pr(y) | ^o)|,=x"-HQ,9o) 

< sup (xVb +lE(Pf (y) |^0)|3^=X"-i(Q,a;)) 

= sup (xyo+E(P['(x"~'(Q,a;))|^o). 

Conversely, setting ggl"^) ■= Iq~^ yields the reverse inequality. 
(5) This item follows from P5. 

5 Affine value function with bang-bang optimal controls 
5.1 The main result 

We recall that, for every integer n > 1, the triangular set of admissible values for a couple of 
(deterministic) global constraint (at time 0), as introduced in P5, is defined as: 

T'^{n) := {{u,v), < u < v < n}. 

Then we will define a triangular tiling of T+(n) as follows: for every couple of integers 

< i < i < n - 1, 

T+. := {{u,v)e [i,i+l]x[jj+l], v > u+j-i} and Tr. ■= {{u,v)e [i, i+l]x[j, j+1], v < u+j-i}. 
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One checks that 

THn) = { (J T;+)U( U ^ij)- 

0<i<j<n-l 0<i<j<n-l 

Theorem 2 The multi-period swing option premium with deterministic global constraints Q : = 
(QmirnQmax)G (u) OS defined by is always obtained as the result of an optimal strategy. 

(a) The value function (premium): 

- the mapping Q i-^ Pq{Q,J^) is a concave, continuous, piecewise affine process, affine on every 
triangle Tr^- of the tiling ofT^{n). Furthermore, 

P^{0, 0) = 0, Po"(0, n) = E ( 1^ V+ I ^0 J and P^{n, n) = E ( |^ | J • 

- If Vi > as. for every i e {0, . . . ,n - 1}, then, for every Q = (Qmin, Qmax) G T+{n), 

Po{Q,T) = Po((0,gmax),^) = PoiiQ maxi ^;max 

a.s.. (in particular Q max ^ PoiiQ max) Vmax 

is a.s. non- decreasing). 

(6) The optimal control: 

- If the global constraint Q = (Qminj Qmax) £ nr+(?i), then there always exists a bang-bang 
optimal control q* = iql)o<k<n~i with is {0, l}-valued for every k = 0, . . . ,n — 1. . 

- IfVi ^0 a.s. for every {0, . . . ,n — 1}, then there always exists a bang-bang optimal control 
which satisfies X]o<fc<n-i = Q max • 

- Otherwise the optimal control is not bang-bang as emphasized by the case n = 2 (see proposi- 
tion\^ below) 

We will first inspect the case of a two period swing contract. It will illustrate in a simpler 
setting the approach developed in the general case. Furthermore, we will obtain a slightly more 
precise result about the optimal controls. 

5.2 The two period option 

We assume n = 2 throughout this section. The first result is the following 

Proposition 1 Let Q = (Qmim Qmax) S T+(2) denote an admis.sible global constraint and Iq := 
[(Qmin — l)"*", Qmax Al] . There is an optimal control q* = iq^^ql) given by 

q*o = argmax.gH {xl^o + lA(Qmax-x)E(y+|^o)-(Qmin-x)+E(yf l^o)} (5.1) 

Ql = lA(Qmax- go)l{yi>0} + (Qmin- go)^l{Vi<0} (5-2) 

so that 

P^iQ,J')=Eiq*oVo + qlVi\To). 

Proof. Let q = iqo,Qi) be an admissible control: qo + qi € [Qmin, Qmax] and qi are [0, l]-valued 
^j-measurable, i = 0,1. Consequently qo is /^-valued and gi is [iQmin — Qo)~^ , (Qmax — 9o)Al]-valued. 
Hence 

90^0 + qiVi < qoVo + lA(Qmax- qo)V^+ - (Qmin" qo^Vf . (5.3) 

On the other hand 

E iqoVo + lA(Qi„ax- ^70)^1+ - (Qmin- qo)^V^~\^o) 

= qoVo + lA(Q„,ax- qoMV+\To) " (Qmin" (Zo)+E(yf |^o) • 
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The mapping a:: i— >■ xVo + lA{Qniax—x)^{Vi~\^o) — {Qmm—x)~^^{Vi \J^o) (called the objective variable 
from now on) is piecewise affine on Iq with ^o-™6asurable coefficients so the above definition of 
defines an ^o-™6asurable /g-valued random variable. Now, combining the above inequalities yields 

E(qo^O + ^1^1 I -^O) < qoVo + lA{Qm..-qoWiVl^\^o)-{Qmin-qo)^^{Vf\To) 

< sup QoVo + lA(Qmax - go)E(y+ I Jb) - (Qmin- go)+E(yf |.7^o) 

= E{q*oVo + qlVi\To). ^ 

In the proposition below we investigate in full details the case n = 2. 

Proposition 2 Let n = 2. The two period swing option premium with admissible global constraints 
Q = (Qmim Qmax) G T+(2) as defined by is always obtained as the result of an optimal strategy. 

(a) The optimal control: 

- If the global constraints Qminj Qmax only take integral values (in {0, 1, 2}) then there always 
exists a {0, l}-valued bang-bang optimal control. When Q simply satisfies Qmax — Qmm £ {0, 1, 2}, 
there always exists a bang-bang optimal control. 

- If Vq, Vi > a.s., then any optimal control q* is a.s. bang-bang and satisfies qQ + q\ = Qmax 
on {Vi > 0, i = 1,2}. 

- Otherwise the optimal control is generally not bang-bang. 

(b) The value function (premium): 

- the mapping Q = {Qmm, Qmax) ^ PoiQi ^) is affine on the four triangles T^j that tile T^(2). 

- Furthermore, when Vq and Vi are a.s. non negative, 

PoiQ,^) = (Qmax-1)+Al {VoAE{Vi\To)) + (Q„,ax A 1) (Fq V E(yi |^o)) • 

The objective variable being piecewise affine on Iq, q^ is equal either to one of its monotony 
breaks or to the endpoints of Iq. Consequently, a careful inspection of all possible situations for 
the global constraints yields the complete set of explicit optimal rules for the optimal exercise of 
the swing option involving the values Vq and ]E(y|^|^o) (expected gain or loss at time 0) at time 
and Vi at time 1. 



Q€ T^Q i.e. Qmin < Qmax < 1: [0; Qmax] and the objective variable reads 

(70 Vb + (Qmax- go)E(F+|^o) " (Qmin " '?o) + E(yf |.Fo) 

with one monotony break at Qmin- One checks that 

ql = Qmax, ql =0 on {Vq > E(y+|^o)}, 

qo ~ Qmin; q\ — (Qmax~ Qmin)l{yi>o} on {E{Vi\J^o) < < ^(Vi'l^o)}, 

qo = 0> 11 = Qmaxl{yi>o} + Qminl{yi<o} on {Vq < E{Vi\J^o)}- 

Note that 

qo — Qmin 5 *?! — Qmax Qmin OH mVi\To) <Vo< E(y+|^o)} n {Vi > 0} 
so that the above optimal control is not bang-bang on this event except if QminS {0, 1} or Qj, 

Qmin- 
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^01 Qmin ^ 1 ^ Qmax ^ 2, Qmin ^ Qmax !• 



Iq = [0, 1] and the objective variable reads 



qoVo + lA(Q^ax- qo)E{V+\To) - {Qmin- qo)+E{Vf {J'o) 
with monotony breaks at Qmim Qmax — 1- One checks that 

q*0 = l, = (Qmax- l)l{Vi>0} on {l/o >E(y+|^0)}, 

(70* = Qmax- 1, gi* = l{y,>o} on {0<yo<E(y+|^o)}, 

'?S=Qmin, ql = l{v,>0} on{-E(yf|^o)<V^O<0} 
Qo = 0, = l{yi>o} + '3minl{\/i<o} on {Vb < -E(yf iJTo)}. 

Note that 

?S + <?!*= Qmax- 1 on {0 < < E(y+ |^o) } H {^i < 0} 

so that the control is not bang-bang on this event, except if Qmax 

G {1,2} or Q 

max — 1 ~l~ Qmin 5 

since the local control (/q and the global constraint are not saturated. Likewise 

q* + ql = l + Q^^ on {-E(yf |.^o) < Vq < 0} R {Vi > 0} 
and the optimal control is not bang-bang on this event, except when Qmin S {0, 1} or Qmax = 

1 ~l~ Qmin- 

Note that both events correspond to prediction errors: V\ has not the predicted sign. Moreover, 
these events are a.s. empty when > a.s., z = 1,2. On all other events the optimal control is 
bang- bang. 



Q€ Tq^ i.e. Qmin < 1 < Qmax < 2, Qmin > Qmax— 1: Then the monotony breaks of the objective 



process (with the same expression as in the former case) still are Qmim Qmax— 1- A careful inspection 
of the four possible cases leads to 

ql = \ = (Qmax- l)l{yi>o} on {Fo >E(y+|^o)}, 

9o — Qmim *?! = (Qmax- Qmin)l{Vi>0} On {E(yi|jP"o) < VJ) < E(T/|^+ | /"q)}, 

ql = Qmax- 1, q{ = l{yi>0} + (Qmin- Qmax + l)l{Vi<0} On {-E(yf jj^g) <Vq< E(Vi|:^o)} 

9o = 0, ql = l{Vi>o} + Qminl{yi<o} on {Vb < -E(yf iJ^o)}. 

Note that on the event 

{-E(yf |.Fo) < ^0 < E(yi|.^o)} n {Vx < 0} 

the optimal control is not bang-bang, except if QmaxS {1, 2} or Qmax = Qmin (both and ql are 
(0, l)-valued) or Qmax = Qmin + 1 (9o = Qmin, ?! = 0); and on the event 

{E(Vi 1.7^0) < Vb < E(y+|.Fo)} n {Vi > 0} 

the optimal control is not bang-bang either (except if Qmin € {0, 1} or Qmax = Qmin or Qmax = 
Qmin + 1) by similar arguments. 

Note that these events do not correspond to an error of prediction. On all other events the 
optimal control is bang-bang. 



Qs T^i i.e. 1 < Qmin < Qmax < 2: The objective variable is defined on Iq = [Qmin — 1, 1] by 



qoVo + lA(Qmax- qoMV+lJ'o) - (Qmin- go)E(yf |.Fo) 
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with only one breakpoint at Qmax— 1- One checks that 

9o = l> = (Qmax- l)l{Vi>0} + (Qmin- 1)1(^1 <0} On {E(T/i | J^q) < ^o}, 

QO = Qmax- 1, ql = l{yi>0} + (Qmin- Qrnax + l)l{Vi<0} On {-EiVflJ^o) < Vq < E{Vi\I^o)} 

qo=Qmin-l, q*i = l on {Fo < -E(yf l-Fo)}. 

Once again on the event 

{-E(yf i^o) < ^0 < ^{Vi\To)} n {Vi < 0} 

the optimal control is not bang-bang, except if Qmax € {1,2} or Q 

max — Qmin 

or Q 

max — Qmin + 1. 

Finally, note that when Vq, Vi > 0, the events on which the optimal controls are not bang-bang 
are empty. ^ 

5.3 Proof of Theorem d 

(a) We proceed by induction on 7i. For n = 1 the result is trivial since T"''(l) = T^^ and Po{Q) = 
<5maxl{vb>0} + Qminl{Vo<o}- (When n = 2 this follows from Proposition 

Now, we pass from n to n + 1. Note that combining the backward programming principle and 
PI yields 

Po"+'(Q,^) =sup{xyo + E(Po"(x"(Q,x),(^i+^W«-i)|^o), xel^}. (5.4) 

We inspect successively all the triangles of the tiling of T^{n + 1) as follows: the upper and 
lower triangles which lie strictly inside the tiling, then the triangles which lie on the boundary of 
the tiling. 



l<i<j<n-l: 



Then, x"(Q,a;) = Q - x(l,l) and = [0,1] . One checks that 

X"(Q,x)GT;+ifxG [0,Qmin-i], x"(Q,a;)G rr^^^. if [Qmin - Qmax-j] and x"(Q,x)G 7;+i^.„i 
if X € [Qmax— j>l] (see Figured]). These three triangles T^j, T~_^- and T^_ij_i are included in 
T"^(n). It follows from the induction assumption that {u,v) i— > Pq{{u,v)^{Tij^)) is a.s. affine 
on them. Hence there exists three triplets of ^i-measurable random variables (A^ , B'^ , C"^) , 
m = 1,2,3, such that, for every Qg T^j, 

3 

^o"(x"(Q,x), (^1+.)) = lj^{A'''{Qmin -x) + B^^\Q^,, -x) + C^) 

m=l 

where Jq = [0, Qmin -i], Jq = [Qmin - Qmax " j] and = [Qmax - j, 1] • Note these random coef- 
ficients satisfy some compatibility constraints to ensure concavity (and continuity). Consequently 

3 

xVo + E(Po"(x"(Q, x), \ro) = Y, lj^{xVo + A^{Q^,, - x) + B^{Q^,, - x) + C^) 

m=l 

where = E{A^\Tq), etc. A piecewise affine function reaches its maximum on a compact 
interval either at its endpoint or at its monotony breakpoints xi = 0, X2 = Qmin — ^5 2:3 = Qmax — J, 
X4 = 1. Hence, 

sup (xyo+E(Po"(x"(Q,^),(^i+.))l-^o)) = max{x,l^o + ^"(Qmin-2;,)+So"(Qmax-x,) + Co"', 

(£,m) = (l,l), (2,1), (3,2), (4,3)}. 
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Figure 1: x i— x"'(^) ^or Qg T^-, l<i<j<n— 1 



It is clear that the right hand side of the previous equahty stands as the maximum of four affine 
functions of Q. One derives that Q ^ Pq^^{Q,!F) is a convex function on T^^ as the maximum of 
affine functions. Hence it is affine since we know that it is also concave. 



Qe Tr., l<i<j <n-l: 



This case can be treated likewise. 



Qe T^j, 1 < i < n - 1: In that case = [0, 1], x'^(Q,x) = ((Q„,in - x)+,Qi„ax -x),xe I^. 

- If g G T+, = Q- x{i,i) £ r+, x e [o,Qmin], x'^iQ.x) = (o,g^ax - x) e t+, 

X £ [(3mm,Qmax - j], x"iQ,x) = (0, Qmax - x) £ To^-i, X £ [Qmax - j, 1] (see Figure E]). The 
induction assumption implies that x i— > Pq'(x"(Q, a;), (^i+.)) is piecewise affine with monotony- 
breaks at Qmin and Q^ax - J- 

- If <5 € Tq~-, x^iQjx) crosses the upper (horizontal) edge of T^j_-^^ at x = Qmax — j and the 
left (vertical) edge of T~^{n) at x = Qmin- Hence x i— > Po{x"iQjx), (.Fi+.)) is again piecewise affine 
with monotony breaks at Qmin and Qmax — j- 

In both cases one concludes as above. 



Qg Tq^o- On proceeds like with Tqj except that Iq = [0,(5max] which yields only one monotony 



break at 



QeT^^^, 1 <i <n-l: Assume first Q G T+^. = [0, 1] and x''{Q,x) = (Qmin - x,n) ii x £ 



0, Qmax - n]. Otherwise x"(Qi x) = (Qmin - x, Qmax -x). It follows (see Figure EI) that x"(Qi x) £ 
i;+ _i if X G [0, Qmin - i] and x'^iQ, x) £ T^t^^^-i if ^ ^ [Qmin - ^, 1]- Both and i;ti„_i are 

included in T'^{n). Hence {u,v) ^ Pq{{u,v), (^i+.)) is affine on both triangles, one derives that 

Po"(x"(Q,x),(-^i+.)) 

= l.„o,Q_-„, (^'(Qmin -X)+ B'n + C^) + (^'(Qmin " x) + B\Q„,^^ - x) + C') 

where A^, S™, C", m = 1,2 are J^i-measurable r.v.. Then, one concludes like in the first case. 
If Qg T-~ , one proceeds likewise except that the two "visited" triangles of T^(n) are T^-^ 



QGTo^: re = [0,1] and 
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- x"(Q,x) = ((Qmin - xe if Qg ro\, 

~ X (Q) {Qm'm -^i J 3^ £ [0, Qniax ^] > X {Qt ■^^ {Qm'm Qmax ? ''^ G [Qr 

n,Qmm] if Tq^. 

In both cases the only "visited" triangle is T^n-i T^{n) and one concludes as usual. 



Q e r+„: = [Qmin - n, 1] and x"(Q, a;) = {Qmm -x,n) ifxe [Q^m - n, Qmax - n], x^iQ, x) 



(Qmin — 2;,Qmax " x) otherwise. Hence x"'{Qjx) takes its values in T^_in_i on which {u,v) ^ 
Pq{{u,v), is affine. The conclusion follows. 

The inspection of all these cases completes the proof of the induction. 

The values of Pq{Q) when Q € {(0, 0), (0, n), (n, n)} are obvious consequences of the degeneracy 
of the global constraints. 

(6) We deal successively with the two announced settings. 

- GLOBAL CONSTRAINTS IN N^: Let n > 1. We rely on the characterization (14. 3|) of 
ql = argmax^g^„-i (xVq + E (Po""^(x"~^(Q, a::), (•^i+^)o<^<n-2) I -^o)) • 



We know from item (a) that {u,v) ^ '^{{u,v),J-i^,) is affine on every tile of T"^(n) 



n—l 



If Q G r+(n) n then one checks that x ^ X^'^iQ, x), x G Iq \ is always affine with Iq 
having and/or 1 as endpoints. To be precise 

-if Q = (ij), 1 < i < i < n-1, x''^HQ,x) = {i-x,j-x)e dT+,^^,ndTr^j^„ I^^' = [0, 1], 

- if Q = (o,i), 1 < i < n - 1, x'^-HQ^x) = (0, j - x) € aro+„i, i^-i = [0, 1], 

,n — „\ fn r\\ 1 



if Q = (0,0), x''~HQ,x) = (0,0), i^-' = {0}, 



ifQ= (i,n), l<i <n-l, = (i-x,n-l)Ga7;ii„_i, = [0,1], 



- if Q = (0,n), x''-HQ,x) = {0,n- 1), /^-^ = [0, 1], 

- if Q = (n,n), = (n - l,n - 1), /^-^ = {1}. 

As a consequence, affinity being stable by composition, x Pq~^{x"'~^{Q,x),Ti+,) is affine 
on/^-i G{[0,1],{0},{1}}. In turn, 

X I— > xVo + ]E(P^~^(x"~"^(Q,x),^i-l_.) I ^o) is affine and reaches its maximum at some endpoint 
of Iq~^ i.e. = or at = 1. Then, inspecting the above cases shows that Q^'* = Q — Qoil, 1) S 
T~^{n — 1) nN^. Using (j4.2p and (14. 3p . one shows by induction on k that is always {0, l}-valued. 

- NON NEGATIVE Vi'. 

Step 1: Global constraint saturated. Let n> 1. Let q* = ('7^)o<fc<n-i be an optimal Q-admissible 
control. We introduce the ^-stopping time 

T{q*) := min {k\qQ-\ \- qI < Qmax - {n - I) + k} 

with the convention min0 = +oo. 

On T{q*) = +00, q^ + ■ ■ ■ + ql > Qmax — {n — 1) + k for every A: = 0, . . . , n — 1. In particular 
the global constraint is saturated at time n — 1, i.e. 

qo + ■ ■ ■ + Qn-l — QraecK- 

Set 

Qk = 9fcl{fc<r(q*)-l} + (Qmax " (n - 1) + T{q*) - (Qq + ' " " + )) l{fc=r(q*)} • 
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One check that g is a Q-admissible control: this follows from the fact that T{q*) is an ^-stopping 
time (note that q* is [0, l]-valued on {T{q*) < +00} since then qQ + • • • + q*-i > Qmax — n + t and 
qQ + ■ ■ ■ + q* < Qmax — {n- "!) + ''")• Likewise one shows that q^ > q^ for every k = 0, . . . ,n — 1. 

Furthermore, note that if Qmax is an integer and q* is {0, l}-valued then q is still {0, l}-valued. 

The Vfc being non negative 

n— 1 n— 1 

J2^kVk>J2qlVk 

k=0 k=0 

hence q is still an optimal control. Furthermore qo ~^ ' ' ' ~\~ Qk — Qmax ~ 

(n- 1) +/c on {/c > T{q*)} 

so that the stopping time r(g) satisfies by construction 

T{q) < T{q*) — 1 on {1 < T{q*) < +00} and T{q) = +00 on {T{q*) = 0} U {T{q*) = +00}. 

Note that if the control q* is bang-bang, iterating the above construction at most n times yields 
an optimal control (7°^* for which T{q"P*) = +00 a.s.. Such a control saturates the global 
constraint. 

As a consequence, this shows that Po{Q,J^) = -Po((0, Qmax), = -Po ((Qmax, Qmax), as- so 
that Qmax I— ^ -Pq ((Qmax, Qmax), -^) is a.s. non-decreasiug and concave. 

Step 2: Local constraints. Since there is an optimal control q* which saturates the global 
constraint, oiiG may assume without loss of generality that ^^min — ^max- 

We proceed again by 

induction on n based on the dynamic programming formula (j5.4p . When n = 1 the result is obvious 
(and true when n = 2 as well). 

Assume now the announced result is true for n > 1. 

Let {0, . . . ,n- 1} and QmaxG + !]• Then, = [0, Qmax A 1] and x'"((Qmax, Qmax), a;) = 

(Qmax ^, Qmax 2;), Iq. HeUCe ^ ((Qmax, Qmax) , 2^) £ '^jj^ ^ [0, Qmax j] and ^ ((Qmax, Qmax), 

Tj'_^j_^, xG [Qmax — j, QmaxAl]. Now V I— > Pq {{v , v) , Ti^,) IS Q.s. concave, non-decreasing, affine 
on [j — and on -|- 1] and non-decreasing. Consequently, there exists i?'", C™", m = 1,2, 
^i-measurable random variables satisfying 

Pb(x"((Qmax,Qmax),x),.Fi+.) = (Qmax - 2;) + C\ X S [0, Qmax - j] , 

-Po(x"((Qmax,Qmax),x),.Fi+.) = (Qmax " 2;) + C'^ , 2; G [Qmax " J, Qmax A 1] . 

with 0<B'^<B^ and B^j + = B^j + a.s.. Set temporarily 

■^{X) := xVo + E (Po(x"((Qmax, Qmax), x),Tl+.) \ ^o) • 

Hence, 

sup ^-(2;) = max(^'(0), ^'(Qmax - i), *(Qmax A 1)). 

Set := E(B'" | ^o) and CJ^ := E(C7™ | Tq) and note that B] < B^ and B^ j + C^ = B^ j + C^ 
a.s.. Elementary computations show that: 

- *(0) < ^(Qmax - j) on {Vo > B^} and ^-(0) > 1'(Qmax - j) on {Vo < B^}, 

- ^'(Qmax - j) < ^(Qmax A 1) On {Vq > B^} and ^'(Qmax - j) > ^(Qmax A 1) On {Vq < B^}. 

Consequently can be chosen {0, Qmax A l}-valued on £^o := {Bq < Vq < -Bq} G and equal 
to Qmax - j on 'Eo := {B^ < Vo < B^} G To. 

On = Eon{q*o = QmaxAljG J-q, one has Po"^n(Qmax, Qmax), ^) = ^'o"^^ ((Qmax, Qmax), 
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Then, the dynamic programming formula shows that the other components {ql.)i<k<n of the 
optimal control on £"0 can be obtained as the optimal control of the pricing problem -Po*(((Qmax — 
l)"*", (Qmax — l)"^)] (-^i+A; H -Eq )o<fc<n.-i ) • One derives from the induction assumption at time n 
that {ql)i<k<n can be chosen bang-bang and ((Qmax — l)"'', (Qmax — l)~'')-admissible which implies 
that q is Q-admissible and bang-bang since q^ = Qmax ^ 1 (011 -^o)- ^ similar proof holds on 

= Eon {q*o = 0}. 

On 'Eo, one has P^'^HiQm..,Qm..),^) = Po^HiQm..,Qra..),^ Eq). Then, the dynamic 
programming formula shows that the other components {ql_^_|^)o<k<n-l of the optimal control on '^Eq 
can be obtained as the optimal control of the pricing problem PQ {{j,j), (^i+fc -E'o)o</c<n-i)- As 
there exists a (j, j)-admissible bang-bang optimal control {ql)o<k<n-i (with respect to 
(^i+A;n'^-Bo)o<fc<n-i) ou '^Eq. Then ql_^_f^ is {0, l}-valued for every fc = 0, . . . , n— 1 (in fact identically 
if J = 0). At this stage one can recursively modify ('?i+fc)o<fc<n-i using the procedure described 
in Step 1 to saturate the upper global constraint. Finally one may assume that X]o<fc<n~i Ik ~ 3 
which in turn implies that {ql)o<k<n is a bang-bang (Qmax, Qmax)-admissible optimal control. ^ 

Application. When a global constraint Q belongs to the interior of a triangle T^, one only needs 
to compute the value of Pq{.,T) at the vertices of this triangle to derive the value of the premium 
at every Q € T^j. When Q is itself an integral valued couple, at most six further points allow to 
compute the premium in a neighborhood of Q. We will use this result extensively when designing 
our quantization based numerical procedure in Section [6l 

An additional result. Proposition [2] shows that it is hopeless to produce in full generality 
some bang-bang optimal control when Qmax — Qram S N. This comes from the fact that at integral 
valued global constraints the bang-bang optimal control may saturate none of the global constraints 
(indeed, so is the case at (2, 2) when n = 2). However, using the same approach as that developed 
in that in the proof of Theorem [21 one can show the following result, whose details of proof are left 
to the reader. 

Corollary 1 Assume the assumptions of Theorem \M hold. If a couple of admissible constraint 
(Qmin,Qmax)G r+(n) Satisfies 

Qmax Qmin ^ {0, . . . , n^ 

then there exist a quasi-bang-bang control in the following sense: W-a.s., q^ is {0, l}-valued except 
for at most one local constraint q^^. 

5.4 The Markov setting 

By Markov setting we simply mean that the payoffs Vk are function of an M'^-valued underlying 
^-Markov structure process {Yk)o<k<n~i i-e- 

Vk = VkiYk), k = 0,...,n-l. 

The Markovian dynamics of Y reads on Borel functions g : M^M 

E{g(Yk+i) I J'k) = E(9(n+i) I Yk) = OkigWk) 

where (0fc)o<fc<n-i is sequence of Borel probability transitions on (R'', ©(M'^)). 

Then the backward dynamic programming principle (j4.ip can be rewritten as follows 

PnQ)=PUQ,Yk), k = 0,...,n 
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with Pn{-, .) = and for every = 0, . . . , re — 1 and every Q G T'^{n—k), 

pUQ,y) = sup{xvk{y) + ekipl+i{x''^-HQ,x),.)){y), xe l^-^-'} , (5.1) 

where iQ,x) = {{Qmin - x)+ , {Q^^, - x)AM) and /^^ - - M)+ Al, Qn^axAl]. 

Pointwise estimation of Pq{Q^). As estabhshed in Theorem [2l one only needs to compute 
the value function Pq{Q) at global constraints Q = (Qmim Qmax) S r+(re) fl i.e. with integral 
components. Moreover, for these constraints, the local optimal control ql is always bang-bang i.e. 
9fcG{0,l}. 

Let = (QjJiijj, Qmax) S T+(re)nN^. For every /c = 0, . . . , re— 1, one defines the set of attainable 
residual global constraints at time k, namely 

QUQ") ■■= {((Qmin - ^) + , (Qmax " ^)+A(re-fe)), £ = 0, . . . , fc} . (5.2) 

(thus Qq{Q^) = {Q^})- Note that the running parameter £ represents the possible values of the 
cumulative purchase process ql. 

One checks that for every Qe Q^iQ^), x"'~^~^{Q, l)x""''"^(<5, 0) € Qfc+i(Q°) since 

X^-'-HQ, 1) = ((Qmin - 1) + , (Qmax " 1) + ) and x"^'~\Q,0) = (Qmin, Qmax A {n-k - 1)). 

Consequently the backward dynamic programming formula having Pq{Q^, y) as a result reads: 

pl{Q,y) = max{xi;fc(y) + efc(p^+i(x"-'=-^(Q, x), .))(y), xG {0, 1} n , Qe QUQ°), k = 0, . . . 

At this stage no numerical computation is possible yet since no space discretization has been 
achieved. This is the aim of the next section where we will approximate the above dynamic 
programming principle by (optimal) quantization of the state process Y. 

6 Computing swing contracts by (optimal) quantization 

6.1 The abstract quantization tree approach 

Abstract quantization In this section, we propose a quantization based approach to compute 
the premium of the swing contracts with firm constraints. Quantization has been originally intro- 
duced and developed in the early 1950' in Signal processing (see [l3]). The starting idea is simply 
to replace every random vector Y : {Q, A) —>■ M.'^ by a random vector Y = g{Y) taking finitely 
many values in a grid (or codebook) T := {yi, . . . ,y7v} (of size N). The grid F is also called an 
iV-quantizer of Y. When the Borel function g satisfies 

\Y-Y\ = d{Y, F) = ^mni^ \Y - y'\, (6.3) 

Y is called a Voronoi quantization of Y (and g as well). One easily checks that g is necessarily a 
nearest neighbor projection on F i.e. satisfies 

Vie {l,...,iV}, {g = y'} c {ueR'^ : \u-y'\= min \u-y^}. 

l<j<N 

The so-called Voronoi cells {g = y*}), 1 < i < -/V, make up a Voronoi tessellation or partition of M"^ 
(induced by F. Note that when the distribution of Y weights no hyperplanes the boundary of 
the Voronoi tessellation of F are Py -negligible so that the Py -weights of the Voronoi cells entirely 
characterize the distribution of Y^ . 
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When p£ [1, oo), the L^-mean error induced by replacing Y by Y, namely 

||y-y|| = 1 E( min ly-yT) 

is called the LP-mean quantization error induced by and its p*^ power is known as the L^- 
distortion. We will see in the next section that the codebook F can be optimized so as to minimize 
the L^-quantization error with respect to the (distribution of) Y. 

Our aim in this section is to design an algorithm based on the quantization of the Markov chain 
(Yfc) at every time k to approximate the premium of the swing contract with firm constraints and 
to provide some a priori error estimates in terms of quantization errors. 

Quantized tree for pricing swing options. As a first step we consider at every time k a grid 
Ffc := {yj., . . . ,y^''} (of size Nk). Then, we design the quantized tree algorithm to price swing 
contracts by simply mimicking the original dynamic programming formula (j4.ip . This means in 
particular that we force in some way the Markov property on (lfc)o<fc<n-i by considering the 
quantized transition operator 

QkigM) = E ■= nn+i = vU, I % = vl) 

so that 

QkioWk) = E {g{Yk+i) I n) , A; = 0, . . . , n - 1. 

[> Let ^ T"^(n)nN^ be a couple of (deterministic) global constraints (at time 0). The quantized 
dynamic programming principle is defined by 

C(Q):=0, QGT+(n) nN^, 

Pfc"(Q) := max(x^;fc(n)+E(P^+i(x"-'=-\g,:E)) | n),x€ a/^-^-^),gG Q^(Q°),fc = 0, . . . ,n - 1.(6.4) 
One easily shows by induction that, for every Qe Q^(Q*') (residual global constraint at time k), 

P^{Q)=K{Q,Yk) 

where KiQ,y) = 0, Q e Q'^.iQ'') , y e R'' , 

K{Q,yl) = sup{x^;fc(yi) + efc(p^+i(x"-'^-^(Q,x),.))(y^),xea/2"i-'=} (6.5) 
i = l,...,Nk, Qe Q^(Q°), k = 0,...,n-l. 

\> When € T~^{n) \ N^, one defines Pq{Q) (and Pq{Q, .)) by affinity on each elementary triangle 
that tiles T+{n). 

Complexity Let us briefly discuss the complexity of this quantized backward dynamic procedure. 
Let k£ {0, . . . ,n — 1}. At every "nod" yl the computation of Qk{Pk+i{x"'~^~^{Q:^): ■))(?/fc) requires 
iV/c_|_i products (up to a constant), so that for a given residual global constraint the complexity at 
time k in the dynamic programming is proportional to Nf^Nf^^i. On the other hand, one checks 
that 

card(Q^(gO)) = (Q^^^Afc) + l - (Q^,, - Q^,, - {n - k) - 1)+. 
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Consequently, the complexity of the computation of p^{Q^ ,Yq) is proportional to 

n— 1 



5^card(Q^(Q°))AffcAffc+i. 



fc=0 

A simple upper-bound is provided by 

n-1 



E((^maxAfc) + l)iVfciVfc+l 



n-1 



k=0 

and a uniform one by 

k=0 

Note that this last upper bound corresponds to the complexity of the quantized version of the 
algorithm based on some penalized global volume constraints (see the companion paper [5]). 

A priori error bounds for the quantized procedure 

Theorem 3 Assume that the Markov process (^fc)o<fc<n-i is Lipschitz Feller in the following sense: 
for every bounded Lipschitz continuous function g : R"^ — > M and every {0, . . . , n — 1}, @k{g) is 
a Lipschitz function satisfying [@k{g)]Lip ^ [©fclLipbjLip- Assume that every function ■ IR'^ — > M 
is Lipschitz continuous with Lipschitz coefficient [vk]up- LetpG [1, oo) such that maxo<fc<n-i \Yk\ E 
LPiJP). Then, there exists a real constant Cp > such that 

n-1 

II sup \P^{Q)-P^mi<CpY.\\^k-Ykl (6.6) 

Remark. In most situations JFg = {0, i}} so that the error term \Pq{Q) — Pq(Q)\ is deterministic. 
When J^Q is not trivial, it is straightforward from (j6.6p (with p = 1) that 

n-1 

sup \E{P^{Q)) - E{P^{Q))\ <CiY, W^k - nlL- 

QGr+(n) fc^O 



We first need a lemma about the Lipschitz regularity of the p^ functions. 

Lemma 6.1 For every k G {0, . . . , n — 1}, the function y i— p^{Q, y) is Lipschitz on , uniformly 
with respect to Q £ T~^{n — k) and its Lipschitz coefficient \p^]up,y '■= supQg'p+(„_;.,)|j»^((5, .)]Lip 
satisfies for every k£ {0, . . . ,n — I}, 

bn-l]Lip,y < K-l]Lip, [Pk]up,y < [vk]up + [S/c]Lip [Pfe+l]Lip,j/. 

Proof. This follows easily by a backward induction on k, based on the dynamic programming 
formula (|5.ip and the elementary inequality | supjgj — supjgj bi\ < supjgj \ai — bi\. ^ 

Proof of Theorem [3l First note that, by piecewise affinity of Pq and Pq, one has 

sup \p^{Q)-p^m= sup \p^{Q)-p^m- 

QGT+in) Q6r+(n)nN2 
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Temporarily set T^{n) := T+(n) n N. Let A; G {0, . . . , n - 1}. Now 

sup \pUQ,Yk) - rk{Q,Yk)\<\vkiYk)-Vkin)\ (6.7) 

QeT+{n-k) 



+ sup 



E(K+i(x"-'-'(Q,x),n+i) i^Ffe) -E(p^+i(x"-i-'(Q,x),n+i) in) 



Now, using that 0^ is a Markov transition and that a{Yk) C iT(yfc)-measurable, one gets 
E(p^+i(x"~'"'(Q,x),n+i)|^fc) - E(i^+i(x"-^-^(g,x),n+i)|n) 

= efc(p^+i(x"-^-nQ, x), .))(n) - E (efe(p^+i(x"-^-*^(Q, x), .)(n) I 

+E x), n+i) - p^+i(x"-'-'(Q, x), n+i) I Yk) 

= efc(p^+i(x""^-^(Q,x),.))(n)-efc(p^+i(x"-^-'=(g,x),.)(n) 
+E (efc(p^+i(x""^"'^(g, x), .){Yk) - efc(p^+i(x"-^"^-(Q, -m) I n) 

+E (p^+i(x"-i"^-(Q, n+i) - i^+i(x""'"'(Q, x), n+i) I n 



Consequently, still using the elementary inequality | sup^g^^^ — sup^,gj5f bx\ < sup^,gj \ax — bx\ 
for any index set X and, for every x€ dl^~^~^, that 

x''-''-HT+{n - k),x) C r+(n - fc - 1) 

(see the proof of Theorem [2)^6)), one has 

sup \pl{Q,Yk)-pk{Q,Yk)\ < \vk{Yk) - Vk{Yk)\ 
QeT+(n-k) 

+ sup |e, (p^+i(Q', .)) (Yk) - Qk{pl+iQ\ .){Yk)\ 

Q'&T+in-k-l) 

+E ( sup \Qu{pl+i{Q\ .){%) - Qk{pl+iiQ', .)(n)l I % 

\Q'6T+(n-fc-l) y 

+E ( sup b^+i(Q',n+i) -p^+i(Q',n+i)l |yj . 

\Q'6T+(n-/c-l) / 

Temporarily set for convenience, A^'^ := || supQg^+(^_^^ bfc(Q,^fc) -Pfe(Q,^A:)lllp- O^^e derives 
that for every k = 0, . . . ,n — 1, 

K'" < (klLip + 2[efc]LipK+i]Lip,,)iin - nil, + Ktv 

Furthermore, A";^-^ < [fn-ilLipH^n-i — ^n-illp- The result follows by induction. ^ 
6.2 Optimal quantization 

Theoretical background In this section, we provide a few basic elements about optimal quanti- 
zation in order to give some error bounds for the premium of the swing option. We refer to [15] for 
more details about theoretical aspects and to [26] for the algorithmic aspects numerical applications. 

Let p€ [l,-|-c»). let Y^ L^(J7,^, P) be an M'^-valued random vector and let > 1 be a given 
grid size. The best L-P-approximation of y by a random vector taking its values in a given grid F 
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of size (at most) N is given by a Voronoi quantizer which induces an L^(P)- mean quantization 
error 



e7v,p(F,r) = ||y-y^||, = Emin|y-y*l^^ 

It has been shown independently by several authors (in various finite and infinite dimensional 
frameworks) that when the grids T runs over all the subsets of M'^ of size at most A^, that eN^p{Y, T) 
reaches a minimum denoted eN^p(Y) (see e.g. [15] or [23]) i.e. the minimization problem 

ejv,p(y) = min|e7v,p(y,r), T C M'*, card(r) < A^} 

has at least a solution temporarily denoted T^^'*\ Several algorithms have been designed to 
compute some optimal or close to optimality quantizers, especially in the quadratic case p = 2. 
They all rely on the stationarity property satisfied by optimal quantizers. In the quadratic case, a 
grid r is stationary 

Y^ = e(y\Y^' 



This follows from some differentiability property of the L^-distortion. For a formula in the general 
case we refer to |16j . In 1-dimension, a regular Newton- Raphson zero search procedure turns out 
to be quite efficient. In higher dimension (at least when d > 3 or 4) only stochastic procedures can 
be implemented like the CLVQ (a stochastic gradient descent, see |23j or the Lloyd I procedure (a 
randomized fixed point procedure, see |13j). For more details and result we refer to |26j . 

As a result of these methods, some optimized grids of the (centered) normal distribution A/'(0; Id) 
are available on line at the URL 

www . quantize .maths-fi . com 

for dimensions d = 1, . . . , 10 and sizes from = 2 up to 5 000. 

It is clear by considering a sequence of grids F^^-* := {r^ , . . . , r^} where {r^)n>i is an everywhere 
dense sequence in that e^^piY) decreases to as A^ ^ oo. 

The rate of convergence of this sequence is ruled by the so-called Zador Theorem (see [30] for 
a first statement of the result, until the first rigorous proof in |15j). 

Theorem 4 (Zador, see U^) (a) Let Ye LP+''(P), p > 1, ?? > 0, such that Py(dn) = ip{du)du + 
u{du). Then 



\\m.NdeN,p{Y) = Jp^d ( / (pp+''{u) du 



1^1 



(6) Non asymptotic estimate (see e.g. ]21^ ): Letp > 1, ?? > 0. There exists a real constant Cd,p,-q > 
and an integer Nd.p.-q > 1 such that or any W^-valued random vector Y , for N > Nd^p^-q, 

eiv,p(y) < Cd,p,r^\\Y\\p+y^N~^. 

Rate of convergence of the quantization pricing method Now we are in position to apply 
the above results to provide an error bound for the pricing of swing options by optimal quantization: 
assume there is a real exponent [1, +oo) such that the (d-dimensional) Markov structure process 
{Yk)Q<k<n-i satisfies 

max iFfclG LP+^(P), n > 

0<fc<n-l 
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At each time A; € {0, . . . ,n — 1}, we implement a (quadratic) optimal quantization grid of Yk 
with constant size N . Then the general error bound result ()6.6p combined with Theorem 21^6) says 
that, if TV > iVd,p,^, 



sup \P^iQ)-P^m 

Q&T+{n) 



< c 



n 



where as the complexity of the procedure is bounded by n(n + l)iV^ (up to a constant). 

In fact this error bound turns out to be conservative and several numerical experiments, as 
those presented below, suggest that in fact the true rate (for a fixed number n of purchase instants) 
behaves like 0{N~d). 

Another approach could be to minimize the complexity of the procedure by considering (opti- 
mal) grids with variable sizes satisfying X^^Zq = nN. We refer to [5] for further results in 
that direction. However, numerical experiments were carried out with constant size grids for both 
programming convenience and memory saving. 



6.3 A numerical illustration 

We considered a two factor continuous model for the price of future contracts which leads to the 
following dynamics for the spot price 

St = Fo,t exp (^(71 £ e-"i(*^^)(iH^i + e-'^^^'-^'^dW^ - ^At^ , t € [0, T], 

where and W'^ are two standard Brownian motions with correlation coefficient p and 

A, = ^ (1 - e-2"^*) + ^ (1 - e-2-^*) + ( 1 - e-("i+"^)0 

2ai ^ ^ 2a2 ^ ^ ai + 02 V / 

Then, we consider a (daily) discretization of the Gaussian process log{St/Fo^t) at times ^ 
(r = 1, n = 365). The sequence (log(S't^/Fo^tfc))o<fc<r!.-i is clearly not Markov. However, adding an 
appropriate auxiliary processes, one can build a higher dimensional (homogenous) Markov process 
(^fc)o<fc<n.-i whom log(5ij./Fo^tj.) is a linear combination. This calls upon classical methods coming 
from time series analysis. Then a fast quantization method has been developed to make makes 
possible a parallel implementation of the quantized probability transitions of y = (yfc)o<fc<n-i- 
For further details about this model and the way it can be quantized, we refer to [5]. In [5], the 
optimal quantization method described above is extensively tested from a numerical viewpoint 
(rates of convergence, needed memory, swapping effect, etc). Its performances are compared those 
of the Longstaff-Schwartz approach introduced in [6]. This comparison emphasizes the accuracy 
and the velocity of our approach, even if only one contract is to be computed and the computation 
of the probability transitions is included in the computation time of the quantization method. 
Furthermore, it seems that it needs significantly less memory capacity when implemented on our 
tested model. 

We simply reproduced below a complete graph of the function Q := (Qmin, Qmax) ^ PoiQ) 
when Q runs over the whole set of admissible global constraints T"^(n). The parameters were 
settled at the following values 

n = 30, ai = 0.21, 02 = 5.4, ai = 36%, (J2 = 111%, p = -0.11 

The graph of the premium function Q 1-^ Pq {Q) defined on T"^(n) is depicted in Figure 1. 

Acknowledgement: Wc thank Annc-Laurc Bronstcin for helpful comments. 



25 




Figure 4: The mapping Q i— > Pq{Q) affinely interpolated from integral-valued global constraints. 
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